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We explore the possibility of performing nuclear structure studies using low energy 
neutrino beams. In particular, low energy beta-beams and conventional sources (muon 
decay- at-rest) are considered. We present results on the total charged- current as well as 
flux-averaged cross sections associated to electron (anti)-neutrino scattering on oxygen, 
iron, molybdenum and lead, as typical examples. It is shown that by using neutrinos from 
low energy beta-beams, information on forbidden states, in particular the spin-dipole, 
could be extracted. 

1. Introduction 

Nuclear isospin and spin-isospin excitations play a crucial role in weak processes. A precise 
description of these transitions is necessary to progress on open issues in high energy 
physics and in astrophysics. Searches for processes beyond the Standard Model, such as 
neutrinoless double beta-decay, require a detailed knowledge of the isospin and spin-isospin 
nuclear response. In particular, beta-decay PQ, muon capture [2113] the 2v double-beta 
decay [I] as well as charge-exchange studies [5|6|7] are essential to constrain neutrinoless 
double-beta decay predictions; while neutrino-nucleus has been pointed out recently as 
a new constraint [SJ. Beta decay is one of the key ingredients of the nucleosynthesis of 
heavy elements (see e.g. PTTU] ). The role of neutrino-nucleus interactions is emphasized 
in [11] for neutrino-nucleosynthesis and for the r-process [T2|13|14p!5] . Such processes are 
also an essential ingredient to determine the response of neutrino detectors which can 
serve as core-collapse Supernova observatories. 

The Isobaric Analogue States (IAS) as well as the Gamow- Teller (GT) giant resonances 
represent the most studied cases. Information on the associated nuclear transitions is 
essentially obtained through beta-decay or charge-exchange reactions. A satisfactory de- 
scription of these excitations is now achieved using microscopic approaches: the shell 
model or the random-phase-approximation (RPA) and its variants. The effective inter- 
actions associated with these models have been extensively discussed and some of their 
parameters have been determined using the allowed transitions as a constraint (see e.g. 
|lb'j). However, one is still not able to explain why the measured GT transition amplitudes 
are always larger than the calculated ones. This is usually corrected by using an effective 
axial- vector coupling constant (the "quenching problem" [PT1[T8] ). 
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Little is still known for the states of higher multipolarity, such as the spin-dipole or the 
higher multipoles. Their importance has been underlined in astrophysics (see e.g. [T3] ) 
and is well established for neutrinoless double-beta decay (see e.g. [TTHI2T)] ). Using sum- 
rules, one can relate the spin-dipole strength distribution to the neutron skin thickness 
defined as the difference between the root mean square radii of the proton and neutron 
distributions. This is an useful information to establish the nuclear matter equation of 
state of neutron stars [2TJ. Some information on these states has been extracted thanks 
to various weak and hadronic probes. A systematic study of muon capture is performed 
in [22] and it is shown that no quenching of the forbidden states is required to account for 
the total capture rates. On the other hand, experiments with (p,n) and ( 3 He,t) have been 
realized [23J. Calculations based on RPA overestimate the experimental results in the 
Gamow- Teller region, while they underestimate them in the 30-60 MeV energy range |21j . 
Therefore if quenching of the forbidden states is necessary is an open issue. This question 
could be settled if more experimental information was available. Furthermore the latter 
could be exploited to improve the parametrization of the phenomenological interactions 
used in microscopic approaches. 

Neutrino-nucleus interactions can probe different spin-isospin and isospin states, de- 
pending on the neutrino energy. While neutrinos in the tens of MeV, such as those 
produced by the sun, essentially explore the IAS and GT excited states, those having 
several tens to hundreds MeV, e.g. from core-collapse supernovea or from accelerators, 
also probe the spin-dipole and higher multipoles. These measurements are in most cases 
inclusive. However, different pieces of information can be obtained by using either sources 
where the energy of the neutrinos can be tuned or complementary probes. So far exper- 
iments are scarce, and include: a measurement on deuteron using reactor neutrinos |24j . 
one on iron |25j and several on carbon with conventional sources (muon decay-at-rest or 
decay- in-flight) [26]. Therefore for the numerous applications, one relies on the theoretical 
predictions available in the literature and based on different approaches. The calculations 
generally agree for neutrinos in the ten of MeV where the response is dominated by the 
IAS and the GT transitions. However, significant discrepancies appear at higher energies 
(see e.g. [271128] for carbon and [30T3T1I3 2 33j for lead). The origin of these disagreements is 
not yet clear and can be due to a different choice in the nuclear structure description (e.g. 
the continuum treatment, the choice of the effective forces, or the correlations included). 

Systematic neutrino-nucleus interaction measurements could be an ideal tool to explore 
the weak nuclear response. At present, new experiments on various nuclei are being 
proposed with a new facility using muon decay-at-rest [31] • Another possibility could be 
offered by beta-beams. This is a new method to produce pure and well known electron 
neutrino beams, exploiting the beta-decay of boosted radioactive ions [35] . The idea of 
establishing a low energy beta-beam facility has been first proposed in [36]. Its potential 
has been discussed in [3"6"|37II38] for nuclear structure studies, in [3"6T3"§] for core-collapse 
supernova physics, and in [4"0~f8|41ll42ll43f44j for the study of fundamental interactions. 
(For a review on the low, medium and high energy scenarios as well as their physics 
potential see [15]). An analysis of how neutrino spectral shapes change at low-energy 
beta-beams depending on the detector geometry and different locations within the same 
detector is made in [IB] , while the possibility of exploiting low-energy neutrinos at off-axis 
from a standard beta-beam is proposed in [4"7] . 
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In this paper, we consider low energy beta-beams as a tool to study the isospin and spin- 
isospin nuclear response. This feature has been first underlined in and explored for 
the neutrino- lead reaction in [35]. Here we perform microscopic calculations of charged- 
current neutrino-nucleus interactions. The nuclear amplitudes are obtained using the 
mean field approach including pairing correlations and the random-phase-approximation 
among quasi-particles. We present the total neutrino-nucleus cross sections as well as 
the contribution of the various multipoles and discuss how their importance evolve, as 
a function of neutrino energy. We give the flux-averaged cross sections associated to 
low energy beta-beams and to conventional sources (muon decay- at-rest), and compare 
the latter to previous calculations, when available. Results are given for oxygen, iron, 
molybdenum and lead taken as sample nuclei. These nuclei present particular interest 
for neutrino detectors and supernova observatories, currently used or under study. In 
particular, oxygen is used in water Cherenkov detectors, like Super-Kamiokande, the 
next-generation UNO [IS], MEMPHYS [H] and Hyper-Kamiokande [50]. An iron-based 
detector was considered in [31]; while molybdenum- [51] and lead [52J ones are being 
studied. Note that the contribution of the different multipoles to the neutral current 
neutrino- lead cross section was studied theoretically in Ref. |53j. 

The paper is organized as follows. The theoretical aspects are described in Section II. 
Results on the isospin and spin-isospin states contributing to the cross sections as well as 
on the total and the flux-averaged cross sections are presented in Section III. Conclusions 
are drawn in Section IV. 



2. Calculations 



The total cross section for the charged current neutrino-nucleus scattering process vi{vi)+ A 
X — > I + A Y is given by 



i / d^rfiEi + Ef-Ev-E^mhyjlHefflvtiK 



(1) 



where Ef (Ei) is the energy of the final (initial) nuclear state, E u (k v ) is the incident 

— * 

neutrino energy (momentum) and E\ (ki) is the outgoing lepton energy (momentum). The 
effective single-particle hamiltonian H e fj is derived by carrying out the Foldy-Wouthuysen 
(FW) transformation for the nucleon weak current operator and retaining terms up to 
0(\q\/M) (M is the nucleon mass and q is the momentum transfer). After performing a 
multipole expansion of the weak current, the differential cross section for the transition 
between nuclear states with angular momenta Ji and Jf is given by: 



da{E u ) G 2 cos 2 9 c E l k l 1 
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where G cos 9c is the weak coupling constant. The a CL and a? represent the longitudinal 
and traverse response of the system to the external field. According to [51] these can be 
expressed as follows: 
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where v = k v /E v is a unit vector denning the incoming (anti) neutrino flux ; /3 —ki/Ei. 
A4j, jCj, Jf and J™ a9 are the charge, longitudinal, traverse electric and traverse magnetic 
multipole operators, respectively. The plus or minus signs correspond to neutrino or anti- 
neutrino scattering respectively. The corresponding matrix elements defined for single- 
particle transitions can be found in |54j . 

A correction to d2J) must be introduced to account for the distortion of the outgoing 
lepton wave function due to the Coulomb field of the daughter nucleus. While the Fermi 
function F(Zf, E{) accounts well for it at low neutrino energies, the modified "Effective 
Momentum Approximation" (EMA) offers a good procedure at high neutrino energies |55j . 
Here we follow the prescription already used in previous calculations, where the Fermi 
function is taken until the corresponding corrections give a larger cross section than those 
obtained with the EMA approximation. 

Flux-averaged cross sections are the measured observables. These are obtained by 
folding the energy dependent cross section ([2]) with the neutrino flux 4>{E V ) that depends 
on the specific neutrino source and the target geometry 



dE v o{E v )<t>{E v 



(5) 



The transition matrix elements entering in Eqs. ([MS]) are calculated within microscopic 
approaches. The starting point is a Hartree-Fock (HF) calculation for the ground-state of 
the nucleus, performed in coordinate space by using the Skyrme-type effective interactions. 
The HF procedure determines the mean-field and single-particle (s.p.) occupied levels. 
The unoccupied levels are obtained by diagonalizing the HF mean-field using a harmonic 
oscillator basis. Therefore, the continuum part of the s.p. spectrum is discretized and 
discrete particle-hole (ph) configurations coupled to a given J n are used as a basis, in order 
to cast the RPA equations in the matrix form. This RPA calculation is self-consistent 
since the residual interaction among ph states is derived from the same Skyrme force as 
the one used to produce the mean field. 

When necessary, to go beyond the closed-shell approximation for the ground state, 
pairing correlations are taken into account in the HF+BCS approximation. Constant 
pairing gaps A p and A n for protons and neutrons are introduced. On top of the HF+BCS 
calculation, the QRPA matrix equations can be written with a procedure which parallels 
what was described above, with the two-quasiparticle (2qp) configurations replacing the 
ph ones. We do not present here the details of the QRPA formalism, which can be found 
in the literature. We simply note that the particle-particle matrix elements are here 
renormalized by means of a parameter g pp that has been chosen to be smaller than 1 
(typically 0.7) to avoid the well-known ground-state instabilities. 

For each multipolarity, every eigenstate of the RPA or QRPA equations is characterized 
by its X* and Y* amplitudes and the transition matrix element for a generic operator 
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Figure 1. Total neutrino- 16 cross section as a function of neutrino energy. Results 
obtained with the SGII and Skyrme III forces are shown. 

O(k) is written as 

(Jf\^20(k)\Ji) — ^2(a\^2d(k)\P)(Xl j3 u a vp + YjipVaUp), (6) 

k a, (3 k 

where a and ft label a given ph or 2qp states, u and v are the BCS occupation amplitudes 
(which reduce to 1 and in the HF-RPA case) and (a\ X)jt 0(k)\ft) are single-particle 
matrix elements. 

3. Results and discussion 

The calculations we show are performed using the Skyrme force. We have checked that 
the use of the SGII force does not change the results significantly (Figured]). The value 
of 

9pp —0.7 is used, but the results are stable against changes of this parameter (Figure 
[2]). All states up to J 71 = 5 ± are included. We present results for four nuclei as typical 
examples, namely 16 0, 56 Fe, 100 Mo and 208 Pb. The total (anti-)neutrino cross sections are 
given in Table [1] for the four nuclei considered. 

We have compared our results to previous calculations, when available, either on the 
total cross sections as a function of neutrino energy, or on the flux-averaged cross sections 
associated to electron neutrinos from decay-at-rest muons. Our results are in reasonable 
agreement with those of Refs. [28.33,29] for oxygen, of Ref. [31.33.57.58j for iron and of 
Ref . [3T1I3 3 59J for lead. As far as the comparison with measurements is concerned, the 
only available is the neutrino-iron cross section. The muon decay-at-rest flux averaged 
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Figure 2. Total (anti)neutrino- 16 Fe cross section as a function of neutrino energy. Results 
obtained for two values of the g pp parameter renormalizing the particle-particle interaction 
in the QRPA calculations (see text) are given. 



cross section is 256 ±108 ±43 (10 _42 cm 2 ). Our result, 352 10~ 42 cm 2 , is in agreement with 
the experimental value. 

The contribution of the different multipoles as a function of the neutrino energy are 
shown in Figure [3] for iron, as an example, both for positive and negative parity states. 
Below 40 MeV the GT J 71 = 1 + transitions dominate the cross section. At 60 MeV the 
J 71 = 1~,2~ become important as well. Beyond 80 MeV all states contribute. 

The details on the neutrino flux associated to low-energy beta-beams are given in [37] 
(Fig. 13]) . Following a preliminary feasibility study, we assume, as in [4"Tf4"2"] . that the 
boosted ions are stored in a storage ring which has 1885 m total length and 678 m 
straight sections, and that the detectors are located at 10 m from the storage ring |56j . 
We consider fully efficient cylindrical detectors as in [37J, namely R = 1.5 m (radius) 
and h = 4.5 m (depth) for metals and R = 4.5 m, h — 15 m water Cherenkov detectors 
for oxygen. The neutrino flux corresponding to conventional neutrino sources, i.e. muon 
decay-at-rest (DAR) is given by the well known Michel spectrum. As seen in Figure S] 
the latter has quite similar shape to a low energy beta-beam with 7 = 10. Note that 
in principle since the cross sections approximately grow as the neutrino energy square, 
the flux-averaged cross sections can show differences due to the high energy part of the 
neutrino flux, as discussed in detail for the case of core-collapse supernova neutrino spectra 

Table [2] presents the contribution (in percentage) of the different states to the flux- 
averaged cross sections Eq.([T]l5]). One can see that in all cases the results for 7 = 10 are 
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Figure 3. Contribution of the positive (left) as well as negative (right) parity states to 
the total neutrino-iron cross section, as a function of neutrino energy. 
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Table 1 

Total cross sections for the indicated neutrino nucleus charged-current reactions as a 
function of incoming neutrino energy. The cross sections are given in units of 10 -42 cm 2 , 
exponents are given in parentheses. 



similar to the DAR case. The neutrino-oxygen cross section is dominated by the 1~, 2~. 
Only with increasing 7, the relative contribution of the 2~ decreases in favor of the 
positive parity multipoles. This behaviour is due to the oxygen structure: positive parity 
transitions require at least 2hu excitations. The situation is different for the neutrino- 
iron, neutrino-molybdenum and neutrino-lead cross sections. For these nuclei the 7 = 6 
case is dominated by the + , 1 + transitions, with a relative ratio of about 1/3. When the 
ion boost increases, their contribution diminishes basically in favor of l - , 2~. For 7 = 14 
the 2 + and 3 + contribution become as important as the Isobaric Analogue state. Table [3] 
presents the results for the anti-neutrino scattering where the anti-neutrino fluxes are 
produced by the decay of boosted 6 He ions. While the 1 + dominates the flux-averaged 
cross section at 7 = 6 for the iron nucleus, about 50% of its value is given by the l - , 2 _ 
transitions at ion boost 14. On the contrary, these states dominate at low gamma for the 
lead nucleus, while almost all multipoles contribute at high gamma. This behaviour can 
be well understood in terms of the underlying shell structure of the lead nucleus, since 
the lowest pn-excitations between the filled and unfilled shells are of negative parity. 

4. Conclusions 

We have presented new results on charged-current neutrino-nucleus interactions. The 
theoretical approach used is the random-phase-approximation among quasi-particles. Four 
nuclei of particular interest have been chosen as typical examples, namely oxygen, iron, 
molybdenum, and lead. We have given both the total cross sections as a function of neu- 
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Table 2 

Fraction (in %) of the flux-averaged cross section associated to states of a given multipo- 
larity with respect to the total flux-averaged cross section, i.e. (a) / '((r) tot . Results are 
given for all positive and negative states having total angular momentum J between 
and 5. The first column gives the considered nucleus, the second tells if the results corre- 
spond to low energy beta-beams or to a conventional source (DAR for the decay-at-rest 
of muons). In the former case the neutrino fluxes are those of 18 Ne ions having Lorentz 
ion boost 7 between 6 and 14. 
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Table 3 

Same as Table [2J but for the anti-neutrino nucleus cross sections and the anti-neutrino 
fluxes produced by boosted 6 He ions. 
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Figure 4. Anti-neutrino fluxes from the decay of 6 He ions boosted at 7 = 6 (dot-dashed 
line), 7 = 10 (dotted line) and 7 = 14 (dashed line). The full line presents the Michel 
spectrum for neutrinos from muon decay- at-rest. 



trino energy and the flux-averaged cross sections associated to possible future facilities 
producing low energy neutrino beams. These are either based on conventional sources or 
on low energy beta-beams. A detailed analysis of the states of different multipolarity (al- 
lowed and forbidden) contributing to the cross sections has been performed, as a function 
of the Lorentz ion boosts. The main goal has been to explore if by combining different 
neutrino ion accelerations, information on the various multipoles can be extracted. We 
conclude that by varying the 7 factor, the role of the spin-dipole states becomes as impor- 
tant as the Isobaric Analogue and Gamow- Teller states allowing nuclear structure studies 
of the forbidden states. 
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